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O : ABSTRACT 

A search for the pair production of doubly-charged Higgs bosons in the process pp — > 
H ++ H- — > /i is performed with the D0 Run II detector at the Fermilab 

Tevatron using inclusive di-muon events. These data taken at an energy of y/s = 

X 
<D 



X 



1.96 TeV correspond to an integrated luminosity of 113 pb 1 and were recorded by 
D0 between August 2002 and June 2003. In the absence of a signal, 95% Confidence 



Level mass limits of M(H^ ± ) > 118.6 GeV and M(i/^ ± ) > 98.1 GeV are set on 



the pair production cross section for left-handed and right-handed doubly- charged 
Higgs bosons, assuming 100% branching into muons and hypercharge \Y\ =2. 



1 Theory 

Doubly-charged Higgs bosons appear in left-right symmetric models, in Higgs triplet 
models and in little Higgs models [2,. The dominant decay modes are expected 
to be like-sign lepton pairs, i/ ±=t — > £ ± £ ± . Decay modes with mixed lepton flavor 
are also possible. In most models, the H ±i: coupling to W pairs is suppressed due 
to the requirement that the vacuum expectation value of the neutral member of the 
Higgs multiplet vanishes. This is needed to constrain the p parameter to unity at 
tree level. 

Pairs of doubly-charged Higgs bosons are produced through the Drell-Yan 
process qq — > j*/Z — > H ++ H . Next-to-leading (NLO) order corrections to this 
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cross section have recently been calculated [3] . The pair production cross sections for 
left-handed states in the mass range studied in this analysis are about a factor two 
larger than for the right-handed states due to different coupling to the intermediate 
Z boson. Left-handed and right-handed states are distinguished by their decays 
into left-handed or right-handed leptons. The cross section also depends on the 
hypercharge Y of the H ±i: boson. 

LEP experiments have searched for pair production of doubly-charged 
Higgs Bosons in e + e" scattering. Mass limits of M(i7 I ±± ) > 100.5 GeV and M(H^ ± ) > 
100.1 GeV were obtained by OPAL T|. and a limit of M(H^) > 99.4 GeV by L3 
for decays into muons. All limits in this Letter are given at 95% Confidence Level 
(C.L.) and for a branching of 100%. Similar mass limits were set for decays into 
electrons |U [5] or r-leptons jU El E] • 

The analysis [7j is based on inclusive di-muon data taken between August 
2002 and June 2003. The total integrated luminosity for the accepted di-muon 
triggers is determined to be 113 ± 7 pb _1 . 

2 Event selection 

The event selection is performed in four steps. The first step (selection SI) requires at 
least two muons, each with transverse momenta pt > 15 GeV, where pt is measured 
with respect to the beam axis. In the second selection (S2), isolation criteria based 
on calorimeter and tracking information are applied to reject background mainly 
from muons originating from semi-leptonic b quark decays. The next selection (S3) 
requires that for events with exactly two muons the difference in azimuthal angle A0 
is less than 47r/5. It is applied to reject Z — > events where one of the charges 

is mismeasured and to further reduce semi-leptonic b decays. This also removes any 
remaining small background from cosmic muons. In the final selection (S4) at least 
one pair of muons in the event is required to be of like-sign charge. These pairs are 
considered candidates for H ±i: — > /x ± /x ± decays. 

When the requirement of having at least one pair of like-sign muons is 
applied at the same time as the selection SI, most of the background from Z decays 
is removed, and only 101 like-sign events remain (Table [TJ. Since no isolation is 
imposed at this stage, the most probable background is due to bb production. 

PYTHIA is used to estimate this background by generating inclusive jet 
events with a minimum transverse momentum for the hard interaction of 30 GeV [Hj. 
The inclusive b quark production cross section a h {p\ > 30 GeV) was measured by 
D0 to be (54 ± 20) nb in the rapidity interval \y b \ < 1 at y/s = 1.8 TeV 0. This 



2 



Table 1: Expected number of background events from Monte Carlo, and number of 
data events remaining after each selection cut. 



Selection 


2 muons 


Isolation 


A0 < 47r/5 


(like- 


•sign) 


p T > 15 GeV 










S4 & SI 


S2 


S3 


Z 


(JL+fi- 


0.9 ±0.3 


0.6 ±0.2 


0.3 ±0.1 


bb 




95.1 ±3.3 


4.4 ±0.8 


0.8 ±0.2 


Z -> 


T + T - 


0.6 ±0.3 


< 0.3 


< 0.3 


ft 




0.24 ±0.01 


0.11 ±0.01 


0.11 ±0.01 


zz 




0.06 ±0.01 


0.05 ±0.01 


0.05 ±0.01 


wz 




0.29 ±0.01 


0.27 ±0.01 


0.23 ±0.01 


MC i 


(sum) 


97.2 ±3.3 


5.5 ±0.7 


1.5 ±0.4 


data 




101 


5 
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cross section, extrapolated to the full y h range and to y/s = 1.96 TeV, is used to 
normalize the bb MC sample. 

In Fig. Q di-muon mass and A0 of the like-sign events are compared to the 
PYTHIA bb simulation. Data and Monte Carlo are in good agreement. Out of 101 
like-sign events, 5 remain after applying the isolation requirement, while 16 events 
remain after applying only the A<p requirement. Assuming that all like-sign events 
originate from bb processes, the isolation and A0 selection efficiencies for bb events 
are 5% and 16%, respectively. Using these the background from bb production in 
the final sample is expected to be 0.8 ± 0.2 events. 

Another potential background are Z — > fi + fi~ decays where one of the 
muon charges is misidentified. We have therefore measured the charge misidenti- 
fication probability in data. 0.45 background events are expected due to charge 
misidentification. This is in good agreement with 0.3 ±0.1 events expected from 
Monte Carlo. Three candidates remain in the data after the final selection. 

3 Limit Calculation and Systematic Uncertainties 

The limit calculations are performed using the program MCLIMIT [10J. It provides 
the confidence level for the background hypothesis, CLb, and the confidence level 
for the signal with background hypothesis, CLs+b, taking into account the expected 
mass distribution for the signal and for the background, and the mass resolution |llj . 
The expected signal rate as a function of the Higgs mass is given by the NLO cross 
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Figure 1: (a) Di-muon mass and (b) A0 distributions for the 101 like-sign events 
remaining in data after the selections SI and S4 (points with error bars), compared 
to the PYTHIA bb simulation (histogram). The five events remaining after the 
isolation selection are shown separately. 

section the signal efficiencies, and the measured luminosity. 

The 95% CL. limit is determined from the confidence level of the signal 
CL S = CL s +b/CL b by requiring CL S = 0.05 1 . 

The following sources of systematic uncertainty affecting the normalization 
of the signal are taken into account: The systematic uncertainty on the luminos- 
ity obtained using the D0 luminosity system is estimated to be 6.5%. The total 
uncertainty on the efficiency amounts to 5%. 

The theoretical uncertainty on the NLO if ±± production cross section 
from choice of parton distribution function and variations of the renormalization 
and factorization scales is about 10% 

The statistical uncertainty on the MC background rate is 27% (Tabled}. 
Adding the systematic uncertainty of 25% on the measured bb cross section [Oj yields 
a total uncertainty on the background rate of 50%. 

4 Results 

Taking into account the systematic uncertainties, a lower mass limit of 118.6 GeV is 
obtained for a left-handed and a mass limit of 98.1 GeV for a right-handed doubly- 
charged Higgs boson, assuming 100% branching into muons, hypercharge Y = |2|, 

x By definition the hypothesis of having a signal plus background is excluded at the 95% C.L. if 
CLs+b < 0.05. 
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Figure 2: Confidence level of the signal, CL S = CL s +b/CL b , as a function of 
the mass M(i/ ±± ) of (a) left-handed and (b) right-handed doubly-charged Higgs 
bosons. The mass regions M(f/ L ±± ) < 100.5 GeV and M(i/"^ =± ) < 100.1 GeV are 
excluded by LEP. The impact of systematic uncertainties is included in the limits. 
The dashed curve shows median expected CL S for no signal. 

and Yukawa couplings > 10 -7 . This significantly extends the previous mass 
limit of 100.5 GeV for a left-handed doubly-charged Higgs boson jl]. 
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